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FINAL  SCIENTIFIC  REPORT 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 

GRANT  NUMBER  AF  AFOSR  83-0338 

'The  Electromagnetic  Launcher  as  a  Circuit  Element' 

Principal  Investigator:  Dr.  Manuel  A.  Huerta 

Department  of  Physics 
University  of  Miami 
Coral  Gables,  Florida  33124 

ABSTRACT 

The  circuit  equations  for  the  inductor  driven  rail  launcher  are  integrated 
numerically  using  a  fourth  order  RuDge-Kutta  method.  Their  solution  is  also 
obtained  as  a  perturbation  expansion  for  small  resistance.  The  EMACK 
parameters  are  used  to  study  the  results.  For  high  current  shots  the  full 
numerical  solution  agrees  with  the  first  order  approximation  to  3%.  For  low 
current  shots  they  agree  to  5%.  A  very  simple  zeroth  order  approximation 
agrees  with  the  full  numerical  solution  to  5%  in  the  high  current  case  and  to 
3%  in  the  low  current  case. 


''  i 


* 
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INTRODUCTION 


To  avoid  repetition  the  initial  proposal  for  this  grant  is  attached  as  an 
Appendix  to  this  Final  Report.  The  work  in  the  grant  was  to  study  the  circuit 
equations  for  the  inductor  driven  rail  launcher.  The  basic  equations  are  Eq. 
(15)  of  the  proposal,  which  comes  from  the  circuit  equation  after  converting 
from  time  to  distance  as  the  independent  variable, 

(Lc  +  L'x)v  ^  +  y~  v2  +  2(R0  +  Rp)  v  +  2R'  J  x  dv  =  b'  , 
and  the  for^e  equation 

m  v  ~  =  L'i2/2  .  (1) 

dx 

This  is  an  integro-dif ferent ial  equation  for  the  speed  v  as  a  function  of  the 
distance  x.  It  is  in  a  form  suitable  for  a  regular  perturbation  expansion  for 
small  resistance. 

We  introduce  dimensionless  variables  in  terms  of  the  characteristic 
parameters  of  the  problem  as  follows.  Let  s  be  a  dimensionless  length  defined 

as 

x  =  gs  ,  (2) 

where 

8  L-  (3) 


We 


introduce  a  velocity  scale  v^  by 


mv , 


I  2 

-  Vo 


(4) 


where  i  is  the  initial  current  in  the  inductor.  We  define  a  dimensionless 
speed  u  by 

v  =  vf  u  , 


1 


and  a  dimensionless  current  I  by 

i  =  4-. 

i  „ 


The  inductance  and  resistance  have  a  time  constant  given  by 


t  =  - ^ — 

1  R  +  R 
o  p 


in  '.eras  of  which  we  define  an  expansion  parameter  e 


Let  k  be  the  resistance  ratio 


R  +  R 

o  p 


and  q  be  the  dimensionless  initial  speed 


q  =  qv0/vf 


where  vQ  is  the  initial  speed  of  the  projectile. 


In  terms  of  the  above  variables,  Eq.  (1)  is  rewritten  as 

( 1  +  s )  u  —  +  ^  +  eu  +  ek  f  s  du  =  b  , 
ds  2  Jn 


f (u2)  -  I2  . 

ds 


where  b  is  given  by 


(1  +  q^ ) 

b  =  -  +  eq 

2 


We  integrate  the  full  Eq.  (1)  using  a  fourth  order  Runge-Kutta  method  and 


take  the  result  as  the  standard  for  comparison  with  the  approximation  scheme 


The  lowest  order  approximation  uq  amounts  to  neglecting  the  resistance 
completely.  Then  the  kinetic  energy  gained  by  the  projectile  is  given  by  the 
drop  in  the  magnetic  energy  stored  in  the  system  inductance  Lq+L'x.  Despite 
its  simplicity,  the  speed  of  Eq.  (14)  agrees  rather  well  with  the  Runge-Kutta 
solution.  It  overestimates  the  speed,  but  the  error  is  not  much  greater  than 
other  errors  that  are  inherent  in  determining  what  are  the  correct  parameter 
values  to  use  to  simulate  the  experimental  situation. 

The  error  is  reduced  if  we  calculate  the  first  order  correction  u^  from 
the  equation 

—  t ( 1  +  s)uoUl]  =  -u0  +  q  -  k  |  0  s  du0  .  (15) 

Q 

After  some  straightforward  integrations  we  obtain 


3 


(k  -  l)(uo  -  q)  [k  +  2 ( 1  +  q  ) ( 1  -k  -k(l+s)] 

“1  =  '  +  T7T  2,172 


2 ( 1+q  )  ’ 


uo(s  +  1) 


klu  (1+s)  -  q] 

+ - bm - 


uo(Hs) 


2 ( 1+q  )  ‘ 

with  the  first  order  speed  given  by 


L(r-uQ)  (r+q)J 


V  =  Vj  (u  +  Uj) 


(16) 


Conclusion 

Below  we  present  the  numerical  results  of  the  full  Runge-Kutta 
integration,  spd  R-K,  compared  with  the  first  approximation 
spd  first  =  vru  ,  and  the  second  approximation 
spd  second  =  Vj.(ao+eUj).  The  results  are  presented  in  the  tables  that 
follow  for  the  EMACK  rail  gun*  shots  1-5.  The  parameters  Lq,  L' , 
D.o+Rp,  R*.  m,  ?r,  etc.  are  given  above  each  table. 

Table  I  contains  the  results  for  shot  1.  It  is  the  least  accurate  case 
because  jf  the  low  current  which  leads  to  a  relatively  large  value  of  e.  In 
this  case  spd  first  is  8%  above  spd  R-K  while  spd  second  is  4%  above  spd  R-K. 
For  the  other  shots  the  current  is  larger  and  the  errors  are  reduced.  The 
value  of  spd  second  is  off  from  spd  R-K  by  only  3%.  It  is  remarkable  that  spd 
first,  which  from  ur  and  neglects  the  resistance  entirely  is  off  by  only  5%. 
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EHACK  shot  number  1 
LO  «  0. 500E-05H, 
RO+Rp  «  0.271E-O4oh«, 
Met  *  0 .  1958+OOKg , 
iO  *  0.500E+06A, 
epsilon  «  .0363936 


LI  *  0.585E-06H,  g  * 

ftl  *  0. 450E-O4ohm,  tl 

barrel  length  *  5.00s 

vO  *  O.OOOE+OOe/s,  vf 


O.lBSE^OOc 
2531 . 8Ss/« 


time 

po£>  l  ti  on 

spd  R-l' 

spd  second 

spd  first; 

0.  OOOOOE+OO 

0.090008+00 

0. 000008+00 

0. OOOOOE+OO 

o.oooooe+ob 

-0.5O637E-03 

0.532^58-02 

0.  6321 2E+02 

0. 63210E+02 

0.632528+02 

■0. 101 27E-02 

0.  21324E-01 

0.  1 26 1 5E+03 

0.  12615E+0^ 

0.  126308+03 

0. 151 91E-02 

0.478998-01 

0.  1 88668+03 

0. 18866E+03 

0.  18901E+03 

0. 20255E-02 

0.  849778-0 1 

0.  25058E+03 

0.  250588+03 

0. 25121E+03 

0. 25319E-02 

0.  13245E+00 

0.31 176E+03 

0.311 77E+03 

0. 31276E+03 

0. 30383E-02 

0.  19017E+00 

0.  37206E+03 

0.372O9E+C3 

0. 37353E+03 

-0.35446E-02 

0.257998+00 

0. 431348+03 

0.43141E+03 

0. 43339E+03 

<J0.40510E-02 

0.  33572E+00 

0.  48946E+03 

0.  48960E+03 

0. 49221E+03 

0.  45573E-02 

0. 42315E+00 

0. 54631E+03 

0.546558+03 

0. 54990E+03 

0.50636E-02 

0.  52006E+00 

0.60177E+03 

0.  60217E+O3 

0 . 60636E+03 

0.  55700E-02 

0. 62621E+00 

0. 65575E+03 

0.65637E+03 

0. 66151E+03 

0.  60763E-02 

0.74134E+00 

0.  70815E+03 

0.  70907E+03 

0.7152BE+03 

0.65826E-02 

0. 86518E+00 

0.  75890E+03 

0.  76022E+03 

0. 76762E+03 

0.  70890E-02 

0.  997438+00 

0.  807938+03 

0.  80977E+03 

0.81B47E+03 

0.  75953E-02 

0.  113788+01 

0.  85519E+03 

0. B5769E+03 

0. B6782E+03 

0.  81016E-02 

0.  128608+01 

0.  90064E+03 

0. 90395E+03 

0.91S63E+03 

0. 86080E-02 

0.  14418E+01 

0.  94426E+03 

0. 94854E+03 

0. 961B9E+03 

0.91 143E-02 

0.  1 6047E+01 

0.  98601E+03 

0. 99146E+03 

0. 10066E+04 

ZO.  96206E-02 

0.  17745E+01 

0.  1 02598+04 

0.  10327E+04 

0. 1049BE+04 

"0.  1O127E-01 

0.  1 9509E+01 

0.  1 0639E+04 

0. 10723E+04 

0. 10915E+04 

0.  1O633E-01 

0.21336E+01 

0.  1 100 IE +04 

0.  1 1 103E+04 

0.  1 1316E+04 

.  0.  1 1 140E-01 

0. 23222E+0 1 

0. 1 1345E+04 

0.  1 1 467E+04 

0.  U703E+04 

0.  1 1646E-01 

0.251648+01 

0.  1 1671E+04 

0.  1 1815E+04 

0. 1 2075E+04 

0.  12152E-01 

0.  27161E+01 

0.  119798+04 

0.  121 48E+04 

0. 1 2433E+04 

£0.  12659E-01 

0. 29207E+01 

0.  1 2270E+04 

0.  1 2466E+04 

0. 12777E+04 

0.  13165E-01 

0.313028+01 

0.  125448+04 

0.  12770E+04 

0. 13108E+04 

0-  1 3671E-0 1 

0. 33441 E+0! 

0.  1 2802E+04 

0.  1 3059E+04 

0.  13427E+04 

0.  14178E-01 

0. 35622E+0 1 

0.  1 3044E+04 

0.  13336E+04 

0.  13732E+04 

0.  14685E-01 

0.  378438+0 1 

0.  1 3271E+04 

0.  13599E+04 

0.  14026E+04 

0.  15191E-01 

0.401028-01 

0.  1 34B4E+04 

0.  1 3850E+04 

0. 14308E+04 

>  0.  15698E-01 

0. 42395E+0 1 

0.  i 36B3E+04 

0.  1 4089E+04 

0.  1 4579E+04 

0.  1 6204E-01 

0.  447208+01 

0.  1 3869E+04 

0. 14317E+04 

0. 1 4839E+04 

0.  16711E-01 

0.470768+01 

0.  140428+04 

0. 1 4534E+04 

0. 15089E+04 

0.  1 7217E-01 

0.  494608+01 

0.  142038+04 

0.  1 4740E+04 

0. 15329E+04 

TABLE  1 
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EMACK  shot  number  2 
LO  -  0.500E-05H, 
RO+Rp  *  t . 271E-04ohm, 
mcs  B  0.225E+O0Kg, 
10  *  0.800E+06A, 
epsilon  K  .0245674 


LI  «  0.585E-06H,  g  «  8.55s 

R1  *  0. 450€-04ohs,  tl  «  0. 185E+00* 

barrel  length  =  5.00s 

vO  =  O.OOOE+OOm/s,  vi  =  3771.24s/* 


¥ 

*• 


time 

O. OOOOOE+OO 
0. 33995E-03 
10.  6799 IE -03 
0.  10199E-02 
0.  1 3598E-02 
0.  1699BE-02 
0. 20398E-02 
C.23797E-02 
,0.  27 1 97E-02 
0.  30596E-02 
0.  33995E-02 
0. 37394E-02 
0. 40794E— 02 
0. 44193E-02 
1 0. 47592E-02 
0. 50992E-02 
0. 54391E-02 
0. 57790E-02 
0.611 90E-02 
,  0. 64589E-02 
10. 67988E-02 
0.  71387E-02 
0.  747B7E-0'?' 
0. 78186E-  02 
0.  B1585E-02 
•  0.  84985E-02 
I  0.  88384E-<'’2 
0. 91784E-02 
0. 95185E-C2 
0. 98586E-02 
0. 10199E-01 
0. 10539E-0 1 
0.  1 0B79E-0 1 
0.  1 1219E-01 
0.  1  1559E-01 


position 

0.  OOOOOE+OO 
0.53386E- 02 
0.21332E-01 
0.4792CE-0: 
0. 85046E -0 1 
0.  1325BF  +  0' 
0. J9040E+00 
0.  25B36Ef  00 
0. 33627E+00 
0.42395E+00 
0.521 16E+00 
0. 6276BE+00 
0. 74326E+00 
0. 86763E+00 
0. 10005E+01 
0. 11416E+01 
0. 12907E+01 
0. 1 4474E+01 
0. 161 14E+01 
0. 17824E+01 
0. i 9602E+01 
0. 21 443E+01 
0. 23346E+01 
r> .  25307E+0 1 
0. 27323E+0 1 
0.29391E+01 
0.31E0BE+0: 
0. 3 367 3E -*-01 
0. 3588 IE *01 
0. 30131E+01 
0. 40420E+01 
f.  .42745E+01 
0. 451 05E+01 
0. 47497E+01 
0. 4991 9E +01 


spd  R-t 

o.  0<  »000E+0G 
0. 94184E+02 
0. 1 8802E+03 
0.2S126E+03 
37370E+03 
4  6509E+03 
0. 55523E+03 
0. 64390E+03 
j. 73093E+03 
0. 81612E+03 
0. 89931E+03 
0 . 98036E+03 
0. 1059 IE +04 
0. 1 1355E+04 
0. 1 2094E+04 
O. 1 2807E+04 
0. 1 3494E+04 
0. 14155E+04 
O. 1 4788E+04 
0. 15394E+04 
0. 1 5973E+04 
0. 1 6526E+04 
O. 1 7051E+04 
0. 17551E+04 
O. 1 B025E+04 
0 .  1 34  7  _'t_  +  >4 
0 .  1  3B97E  i  .'4 
0. 1 9297E+04 
0. 19675E+04 
0. 20030E+04 
0.  2‘  '363E-m;>4 
0 . 20676E  +  04 
O. 20970E+04 
0.21 244E+04 
O. 21501 E+04 


spd  second 

O. 00000E+00 
0.941 S3E+02 
0.  1 8802E+03 
0. 28 1 27E+03 
O.  37370E+03 
0.  4651 IE +03 
0. 55528E+03 
0.  64401E+07 
O.  731 13E+03 
0.  81647E+03 
O.  89989E+03 
0. 9B125E+03 
0.  1 0605E+04 
0.  1 1374E+04 
0.  12121E+04 
0.  1 2B44E+04 
0.  1 3543E+04 
0.  14217E+04 
0.  14868E+04 
0.  15494E+04 
0.  16096E+04 
0.  16675E+04 
0.  1 7230E+04 
0.  1 7763E+04 
0.  1 G274E+04 
0.  1 S763E+04 
0.  19231 E+04 
0.  196782+04 
0.  201 07 E  +  04 
0.  20516E+04 
0.  20907E+04 
0.2 1201 E+04 
0.  2163SE+04 
0.  21 979E+04 
0.  22304E+04 


spd  first 
0. OOOOOE+OO 
0. 94224E+02 
0. 1BB17E+03 
0. 28162E+03 
0. 37433E+03 
0.46610E+03 
0. 55671E+03 
0. 64599E+03 
0. 73374E+03 
0.81982E+03 
0.90409E+03 
0. 98641 E+03 
0. 10667E+04 
0. 11449E+04 
0. 12208E+04 
0. 1 2945E+04 
0. 13660E+04 
0. 14351 E+04 
O. 15020E+04 
0. 15666E+04 
0. 162B9E+04 
0. 168B9E+04 
0. 1 7468E+04 
0. 18025E+04 
0. 1 8561E+04 
0. 1 9077E+04 
0. 19572E+04 
0. 20049E+04 
0. 20507E+04 
0. 20946E+04 
O. 21 369E+04 
0. 21775E+04 
O. 221 65E+04 
0. 22540E+04 
0. 22899E+04 


TABLE  2 


EWACK  \  ■  »  •  number  7 
LO  «  0,i~OE-05H, 
RO+Rp  «=  0.27lE-04ohm, 
Mks  *  0.275E+00Kg, 
iO  =  0.120E+07A, 
epsilon  =  .0181069 


LI  =  0.585E-06H,  q  « 

Rl  =  0. 450E-04ohm,  tl  = 

barrel  length  «=  5.00m 

vO  =•  0.  OOOE+OOm/fc,  v-f  = 


8.  55m 
0. 185E+00S 

51 16. 82m /s 


time 

).  OOOOOE+OO 
).  25056E-03 
:i.  son  ie -or 

:>■  75166E-07 
1 0022E-02 
5.  1 2528E-02 
j.  15034E-02 
0.  17539E-02 
3.  20O45E-02 
0. 22550E-02 
0.25055E-02 
0. 27561E-02 
0. 30066E-02 
0. 32571E-02 
0. 35077E-02 
0. 37582E-02 
0. 4008BE-02 
0. 42593E-02 
0. 45098E-02 
0. 47604E-02 
0.50109E-02 
0. 52615E-02 
0. 551 20E-02 
0. 57625E-02 
0. 60 1 31 E— 02 
0. 62636E-02 
0. 651 4 IE— 02 
0. 67647E-02 
0. 701 54E-02 
0. 72660E-02 
0. 75 1 67E-02 
0. 77674E-02 
0. 801 80E-02 
0. 826B7E-02 


position 

•'i .  O0O00E +  00 
0.533922-  02 
•I  .  2133T-  0  3 
O.&JVSZZ-'.  3 
i  _  S50E3E  — *  : 
0 .  1 32 65L  + 

0. 3  9053E+00 
0. 25956E+00 
t.  .  33657E  +  00 
0. 42438E+00 
0. 52175E+00 
0. 62847E+00 
O. 74429E+00 
0. 86895E+00 
0. 10022E+01 
0. 1 1437E+01 
0. 1 2972E+0 1 
0. 3  4504L +01 
0. 1 6150E+01 
0. 178672+01 
0. 1 9652E+01 
0. 21501E+01 
0. 234 13E+01 
0. 25384E  +  0 1 
0. 2743 0E+O3 
0 . 2'v4r  OE  +  0  3 
0.  33  623  2+91 
0. 3379°E+03 
0. 36022E+01 
0. 332S7E+01 
0. 40597E+01 
0. 429302 i Ol 
0. 453152+01 
0. 47727E+0 1 


spd  R-f. 

;  OOOOE+OO 
O. 1 276 IE +03 
0. 255 1 BE +  03 
0. 3B1S1E+03 
•1  .  E0~3o2+  03 
0. u3 1 55E+ 03 
0.  *,5409E  +03 
O. 874682+03 
C, 993O8E+03 
0. 1 1090E+04 
0. 12223E+04 
O. 1 3328E+04 
0. 14402E+04 
0. 1 5444E+04 
0. 1 6453E+04 
0. 1 7 42BE+04 
0.  1836T+<:4 
0. 1927 IE +04 
O.  23+ *  1  39E+  04 
0. 2097CE+04 
0. 21766E+04 
0. 22525E+04 
O. 23243E+04 
0.  2397-72+  04 

24  5932+04 
0.2521 12+ 04 
0. 257982+ 04 
0.267532 i 04 
0. 26877E+04 
0.273712+04 
*0.  2' '’8702+  3  4 
0.252742*04 
O.  236352 +04 
0 . 290702+  *.’4 


spd  second 

0. OOOOOE+OO 
0. 127B0E+03 
0.2551 BE +03 
0. 38  1  SlE+03 
0. 507372+03 
0. 6315BE+03 
0. 754 15E+03 
0. 874B2E+03 
0.99335E+03 
0. 11095E+04 
0. 12231E+04 
0. 1 3340E+04 
0  1 4420E+04 
0. 15470E+04 
0. 164B9E+04 
0. 17476E+04 
0.  3  8432E+04 
0. 1 9355E  +  04 
0. 20246E+04 
0.21 104E+04 
0. 21 930E+04 
0. 22725E+04 
0.23489E+04 
0. 242222+04 
0.24926E+04 
0. 256012+04 
0. 26247E+04 
0. 26C672+  04 
0. 27460E+04 
0. 2B028E+04 
0. 2B57 1 2+04 
0. 29091E+04 
0. 29539E+04 
0. 30065E+04 


spd  first 
O. OOOOOE+OO 
0. 12784E+03 
0.255342+03 
0. 3B216E+03 
0. 50800E+03 
0. 63257E+03 
0. 75559E+03 
0.87680E+03 
O. 99597E+03 
0. 1 1 129E+04 
O. 12273E+04 
0. 13392E+04 
O. 14482E+04 
O. 15544E+04 
O. 1 6576E+04 
0. 1 7578E+04 
O. 1 8549E+04 
0. 1 94B9E+04 
O. 20398E+04 
0. 21276E+04 
0. 221 23E+04 
0. 22940E+04 
0. 23727E+04 
0. 24485E+04 
0.252142+04 
0. 25916E+04 
0. 26590E+04 
0. 27239E+04 
0. 27B62E+04 
0. 28460E+04 
O. 29036E+04 
0. 29588E+04 
0. 301 19E+04 
0. 30629E+04 


TABLE  3 


EMACK  shot 

number  4 

’  "1 

LO  «=  0.500E-05H,  LI 

=  0. 585E-06H, 

g  =  8 . 35m 

RO+fip  =  C  . 

271E-04ohm,  R1 

=  0.450E-04ohm,  tl  =  0.  lL5E+00s 

■ass  *  0.317E+00Kg,  barrel  length  = 

5.  00m 

iO  «=  0.210E+07A,  vO 

=  0.000E+00®/ 

s,  vf  -  8 

740. 17m/s 

|  epsilon  = 

.0111088 

: 

time 

posi  t  i  on 

spd  R-l 

spd  second 

spd  first 

O.OOOOOE +0'' 

0. OOOOOE 4  00 

0 .  O0000E  4  no 

O. .  i  >OOOOE+O0 

0. OOOOOE+OO 

0. 15372E-03 

O. 53798E-02 

0.  2'.»836E407 

V . 20336E+0 3 

0. 20840E+03 

0. 3074  4E-07 

0. 21342E-01 

0.  4  1 608E+O3 

0.41 608E+03 

0.  41624E+07 

; 

[  0.46116E-07 

0. 47960E- 0! 

7.622fc5E407 

0. 622o6E+07- 

0.  62301 E +<03 

o.  <si4eeE-  " 

0.B5127E-V2 

•>.  22 75784  33 

: .S275BE+03 

..  .  G2B21E+07 

C.  7686 IE-  ' 

.  13Z73EJ  00 

0.  1  .>7'.»7E4  :  4 

n .  1  ;:704E  +04 

0. 10714E+04 

0.92234E-C7 

0.  1 90  67  E  +  O'.1 

0. :2305E404 

0. 127068+04 

0.  12320E+04 

0, 10761E-02 

0.  258780  00 

0.  142758-m'4 

0. 14277E+04 

0.  14297E+04 

_  0. 12298E-02 

0. 7369OE+0O 

0.  1  till  1  E+o4 

0. 16715E+04 

0.  168 41 E+04 

A  0. 13335E-02 

0. 424S4E+00 

0. ISlOOE+O-l 

0. 38. 15E+04 

0. 1B149E+04 

i 

0. 15372E-02 

0. 52279E+00 

0.  19962E+04 

0.  1  9^  74E+04 

0.  20016E+04 

i 

0. 16909E-02 

0. 62933E+00 

0. 21771E+04 

0. 21790E+04 

0.  2 1841 E+04 

i 

0. 18446E-02 

0.  74542E+00 

0.  23531 E+04 

0. 23559E+04 

0.  22  4 22E+04 

\ 

0. 199S3E-02 

0. 87039E+00 

0.  25239E+04 

0.25281 E+04 

0. 25355E+04 

0. 21 520E-02 

0. 10040E+0 1 

0.  26E95E+04 

0. 26953E+04 

0.  27040E+04 

- 

0.23057E-02 

0. 1 1459E+01 

0.  29496E-KI4 

0. 2B574E+04 

0.28676E+04 

0.24594E-02 

0. 1 2959E+0 1 

0. 70041 E+04 

0. 30144E+04 

0.  30262E+04 

* 

0. 261 31E-02 

0. 14537E+01 

0. 315268+04 

0. 31663E+04 

0. 31 798E+04 

0. 27669E-02 

0. 16189E+01 

0. 32957E+04 

0. 331 27E+04 

0.  332B2E+04 

0. 29206E-02 

0. 17913E+01 

0. 34329E+04 

0. 34544E+04 

0.  34717E+04 

j  0. 30743E-02 

0. 19706E+01 

0. 35642E+04 

0. 35907E+04 

0.  36101E+04 

■  0.322B0E-02 

0.21565E-01 

0 . 36B97E+04 

0. 37220E+04 

0. 37436E+04 

0. 33817E-02 

0. 274 86C < 01 

C-.3C5O95E+04 

0. 784838  +  04 

0.  38722E+04 

■  0. 35354E- 02 

0.  25468L+  01 

0 . 39236L'4  04 

0. 7967 78 4  04 

0.  39960E+  04 

0. 36891E-02 

0.  27506E+C  1 

0.  40722E+04 

u.  40S63E+O4 

0.41 153E+04 

i 

0. 3B428E-02 

0. 295?r'84  •; 

.  41353E-04 

u  •  4  1  43 .  i_T  '.>4 

0. 42299E+04 

a  0. 39965E-02 

0.  3 17431 J  01 

.  42331  £  >-0-1 

0. 470578+04 

0.  43402E+04 

9  0. 41503E-02 

0.73976E  +0  3 

•:  .  4  72688-!  0  4 

0. 440898+04 

0.  44462E+04 

0. 43040E-C2 

0 .  36 1  /6E +  O 3. 

i.  .  4-'  1  04L+04 

0. 450778+04 

0.  45481 E+04 

0.44578E-02 

0. 384588+0 1 

0. 447638+04 

0. 460268+04 

0. 46461E+04 

0.461 16E-02 

0.  4078700 1 

0 . 4-'67448'  04 

0. 46°34E+04 

0.  47401E+O4 

-  * 

0. 47654E-02 

0. 431468 i 01 

0.  *1 6480E4  4  4 

0. 47306E+04 

0.  4B305E+04 

•  v 

0. 4r 3  92E-02 

C.45545E+01 

'I  .  -717  1L+04 

0. 48640E+04 

0.  49 1 74E+04 

0. 50730E-02 

0. 479798 +  01 

0. 4"EC6C+04 

0. 49441 E+04 

0. 5000EE+04 

TABLE  4 
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EMACJ  shot  nuir.to'  5 
LO  «=  0.500E-05H, 

RO+Rp  *  0.271E-04ohm, 
mbs  *  0. 3l8E+00Kg, 
iO  =  0. 160E+07A, 
epsilon  «=  .0146033 


LI  -  .  585E-06H,  Q  - 

R1  -  0. 450E-04ohm,  tl  = 
barrel  length  =  5.00m 

vC  -  0. OOOE+OOm/s,  vf  ; 


8.  55m 
0.  185E+00s 

6344 . 41m/s 


time 

P DEI 1 1  on 

cpd  R-l 

spd  second 

spd  first 

0. 00000E-*  00 

(  1  _  1  1M.  1  |<  '2  4  •_»<_■ 

"  ’  >OoOE -*  00 

O.OOOOOE+OO 

O.OOOOOE+OO 

0, 20208E-03 

0. 533"'5C-02 

.  ■.  3  5S49E+*.  7 

0.  158492+07 

0.  15B53E+03 

C.40415E-03 

<  2] 7""E-0: 

•:  .  33  6462  *03 

0.  31 646E+03 

0.  31662E+03 

0.60622E-07 

0.4TP52E  c: 

O.m-353:  -l 03 

0.  47353E  +  03 

0. 473BBE+03 

0. B0B30E- 07 

■•.£•51  032  -0 1 

:  . 1 29312 *  03 

0.62932E- 03 

0.  62975E  +  03 

0. 101 04 E- 02 

0.  3  326 9E  4  00 

0 .  ~3742E+03 

0.  7C34faE+03 

0. 78445E+03 

0.  123  25E- 02 

‘  .  1  9060E+0'. 

0. 9325 1E+03 

0.  93559E+03 

O.  93703E+03 

0. 1 4 1 45E-02 

i.  .  21S£  ’’E 

0. 108522+04 

0. 1 0054 E +04 

0.  10874E+0- 

0.16;  66E-01 
C.  131 B7E-C2 
0. 20207E-C2 
0. 22228E-02 
0. 24249E-02 
0. 26269E-02 
0. 2B290E-02 
0. 3031 0E-02 
0.3277 IE-  02 
0.34352E-02 
0. 36372E-02 
0. 38393E-02 
0.  404 13E-02 
0. 42434E-02 
C.  44455E-02 
0.  464"'5E-02 
0.  48496E-02 
U . 505 1 6E-02 


..  .  33t74F  +  '.X' 
0.  4 24* 12+ 00 
■j .  52207E-'  O'. 
C . 62870E+CO 
0. 744B6E+00 
0. 06967E+00 
0.  1003  IE -*01 
0. 1 1443E+01 
0.  1 29452+01 
0. 1 4520E+0 1 
0. 163  £?E  *  01 
0 .  1 7870E+  01 
0.  19679E  *-01 
0. 21533E+01 
C. 27450E+01 
0 . 25476E+01 
.  ?>7r_?E+01 
>.  2734-3E+01 


0. 

0. 

0. 

c. 

0. 


1  2723E-*  04 
137633 ‘ 04 
15171E+04 
15543E+  04 
17B7BE+04 
191 75E+04 


0. 64666E— 02 
0.66627E  02 


0. 20470E-*  04 
0. 21643E+04 
0. 228172+04 
0.  23?3r,2  *-04 
0. 2502 1  E-*  04 
0.  260520  04 
0.  27050E+  04 
0. 27990C+04 
0.  289022  -04 
0.29763E '04 
0.  3028  1E+04 
0. 313582+04 


0.  12326C+  04 
0.  1 3769E+04 
0.  151B0E+04 
0.  1 6550E+04 
0.  1 7900E+04 
0.  19206E+04 
0.  20474E+04 
0.  2 1 703E+04 
O. 22892E+04 
0. 24042E+04 
0. 25152E+04 
0. 26222E+04 
0. 27253E+04 
0. 28245E+04 
0. 29199E+04 
0.  301 1  52  *-04 
0. 30995E+04 
0.  31S77Z+  04 


0. 1 2752E+04 
0. 1 3B02E+04 
0. 1 5222E+04 
0.  16610E-*04 
O. 1 7963E+04 
0. 192B1E+04 
0. 20562E+04 
0. 21805E+04 
0. 23010E+04 
0. 24177E+04 
O. 25305E+04 
0.26395E+04 
0. 27447E+04 
O. 28461E+04 
0. 29438E+04 
0 . 30379E-*  04 
0. 312B4E+04 
O . 32 1 55E+04 


|  0. 52537E-02 

C  .  316222-’  01 

C. 320° 'E 1 04 

0. 326492+04 

0.329732+04 

0.54558E-02 

0.333672+01 

0.327912  >-04 

0.  3342S2-l  04 

0. 33798E+04 

0. 565B0E-02 

0. 36099E+01 

0. 334492+04 

0. 34 16 ^2+04 

0. 34572E+04 

0.5860 IE-02 

0. 38373E+0 1 

0 . 34070E  *04 

0. 34832E+04 

0.  353 1 6E+04 

0. 60623E-02 

0.  4G£  822+01 

•'  .  346562+04 

0. 318642 '  CO 

0.  36030E+04 

0.62644E-02 

0. 4304 12 • Cl 

0  .'35207E  •  04 

0.  362 152  +  '.-4 

0 . 367 1 6E+04 

).  4343‘‘'E+0: 
*.47C53Ft01 


35~r25E+04 

30213C-!C^ 


0. 363442- 
0.374  141- 


04  0.37376E+04 

04  0.3B010E+04 


TABLE  5 
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Introdoot ion 

This  proposal  is  •  follow-up  to  the  rssssroh  tbs  P.I.  performed  at  the 
AFATL/DLDL  at  E|lin  AFB  under  tha  1982  USAF-SCEEE  Sumer  Faculty  Research 
Program.  That  research  identified  the  basic  research  iaaues  involved  in 
improving  the  models  for  the  arc  plasma  produced  in  electromagnetic  rail 
launchers  with  plasma  driven  projectiles.  The  results  of  that  research  are 
summarized  in  the  final  report  entitled  Basic  Research  Issues  in 
Electromagnetic  Rs il  Lannchers  with  Plasms  Driven  Pro i ec t i les .  A  copy  of  that 
report  is  attached  as  the  Appendix  to  this  proposal. 

The  small  scale  of  s  minigrant  dictates  that  a  problem  of  the  appropriate 
scope  be  considered.  I  propose  *.o  study  the  behavior  of  rail  launchers  as 
circnit  elements.  *iie  circuit  propertied  most  directly  describe  how  the 
stored  electrical  energy  is  converted  into  kinetic  energy  of  the  projectile 
and  determine  the  overall  efficiency. 
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Previous  fork  fej  QlMt  A&tJfeP-H 

The  current  variations  in  rail  launohers  ara  predominantly  of  lcnr 
frequencies.  Thia  anablaa  oa  to  limit  oar  traataant  to  tha  caaa  where  tha 
length  of  tha  raila  ia  anch  ahortar  than  tha  wavelength  of  tha  electroaagnetic 
fielda  preaant.  Under  these  conditions  the  current  is  nnifora  as  a  function 
of  position  along  tha  rails  and  the  distributed  capacitance  is  negligible.  Ve 
ara  than  able  to  speak  of  tha  relationship  between  tha  voltage  at  the  breech 
and  tha  currant  through  tha  rails  and  plasna. 

A  further  s iaipl  if  icat ion  is  possible  regarding  tha  skin  depth.  It  is 
shown  in  the  Appendix  that  provided  that  the  rail  thickness  is  much  smaller 
than  the  rail  separation  one  may  take  the  rail  resistance  per  unit  length  R' 
and  the  inductance  per  unit  length  L'  to  be  constants  independent  of  the 
frequency.  From  Eq.  (36)  of  the  Appendix  we  get  that  the  breech  voltage 
Vb(t)  is  related  to  the  current  i(t)  by 

Vb(t)  -  i(t)  Rb  +  d/dt(Lt,i)  (1) 

whera  Rb  *=  Rp  +  R'x  and  Lb  **  L*x,  (2) 

with  Rp  equal  to  the  plasma  resistance.  Here  x  is  the  distance  from  the 
breech  to  the  plasma  as  shown  in  Fig.  1  of  the  Appendix.  The  i-V  relation 
involves  the  variable  x.  The  equation  for  it  is 

m  d2x/dt2  -  1/2  L* i2  -  f.  (3) 

1  2 

This  equation  is  the  statement  that  the  magnetic  force  j  L'i  accelerates 
the  mass  m  against  a  retarding  friction  f. 

The  form  of  Eq.  (1)  can  be  understood  siaq>ly.  The  voltage  drop  at  the 
breech  is  made  up  of  the  resistive  drop  plus  the  rate  of  change  of  the 
magnetic  flux  #  -  L^i  »  L'xi.  The  fact  that  the  force  on  the  armature  is 
1/2  L* i2  can  be  understood  by  considering  the  energy  into  the  breech.  The 


p07«r  into  tli*  br**oh  It  given  by  iVb .  Fro*  Eq.  (1)  on*  readily  get*  that 

iVb  -  i2!*  ♦  d/dt  <1/2  L'xi2)  ♦  1/2  L'i2  dx/dt.  (4) 

¥•  •••  that  the  po**r  into  tb*  br**oh  go*t  into  *  resistive  h**t  lots,  an 
increase  of  the  aagnetic  energy  (1/2  L^i  )  stored  in  the  rails,  and  a 
third  term  that  gives  the  rat*  of  doing  work  on  the  armator*.  From  this  tern 
«•  deduce  that  the  magnetic  force  on  the  armator*  it  given  by  1/2  L’i  . 

The  circuit  behavior  of  the  rail  launcher  is  described  by  Eqs.  (l)-(3). 
Of  course  the  rails  will  be  part  of  a  coaplete  circuit.  The  simplest  launcher 
arrangenent  consists  of  a  storage  inductor  L0,  with  its  resistance  R0. 
connected  in  series  with  the  breech  as  shown  in  Fig.  1.  Previous  steps  have 
energized  the  inductor  L0  so  that  at  t  ■  0  the  current  has  the  value  i0. 


Fig.  1  Schematic  of  rail  gun.  Vb  is  the  breech  voltage.  The  arc  has  moved 
to  a  distance  x  fron  the  breech, 
resistance  Rc. 


L0  is  a  storage  inductor  vith  internal 
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The  complete  oircoit  equations  are  than 

d/dt [(L0  ♦  L'x) il  ♦  (B*  ♦  Rp  ♦  l*x) i  -  0, 

and  d2x/dt2  -  L'i2/2n. 

where  the  foroa  of  friotion  haa  bean  neglected  in  eq.  (6). 
conditions  are 


(5) 

(6) 

The  initial 


i(t  -  0)  -  i0.  x(t  -  0)  -  0, 

and  dx/dt(t  ■  0)  •  v0  .  (7) 

The  initial  velocity  allova  for  a  preaccelerated  projectile  that  is  injected 
into  the  launcher  with  a  speed  v0.  The  aysten  of  Eqs.  (5)  and  (6)  is 
nonlinear  and  has  not  been  solved  analytically  in  nontrivial  cases.  Nuaerical 
resnlts  are  easy  to  obtain.  In  fact,  Mr.  Kenneth  Cobb,  of  AFATL/DLDL,  already 
haa  a  program  to  integrate  then. 


Previous  wort  by  the  LL 


Scat  of  the  work  under  this  proposal  wonld  50  into  developing  further  the 
nuaerioal  results.  The  sisin  effort,  however,  will  go  into  developing 
analytical  results.  Analytical  results  are  essential  to  understand  how  the 
solution  scalos  as  the  various  paraaeters  are  varied.  The  P.7,  has  already 
a>ade  sufficient  progress  on  this  probleu  that  valuable  results  are  aluost 
assured. 


Ve  first  note  that  uultiplying  Eq .  (5)  by  i  and  using  Eq.  (6)  yields  a 
single  third  order  equation  for  x(t). 


d3x 


dx  d2x 


d2x 


(L°  +  L'x)  J?  +  2L'Tt  d?  +  2(R°  +  +  R'x>  dF  “ 


0. 


(8) 


This  equation  is  to  be  solved  with  the  boondary  conditions  of  Eq.  (7)  together 
with 


0) 


(9) 


Ve  note  that  Eq.  (8)  can  be  rewritten  as 

S  K  +  £  ]  *  r  J 7  *  V  lTt  [57]  ♦  Sr  J 


dt" 


dt 


Obviously  Eq.  < 1 J )  can  be  integrated  once  if  the  R'  tera  is  negligible.  This 
tera  would  be  zero  in  an  ideal  rail  launcher  with  supercooled  rails.  Ve  will 
carry  along  the  R'  integral  tera  and  treat  it  as  a  perturbation.  Integrating 
Eq.  (10)  yields 
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<L, 


d^x  L'  rdxi1  dx  f 

*  L#a)  i?  *  r  y  *  2<tt°  +  V  ST  * 2R'  [ 


d2x 

1? 


dt  - 


♦  2(Ro  «•  Rp)v0 


■  b' 


(11) 


Eq.  (11)  follows  immediately  f row  conservation  of  energy,  that  is 

\  L0i02  f  \  mv02  -  j  (L0  +  L'x)i2  ♦  \  mv2  *  J  i2(Ro  +  Rp  ♦  Rix)dt.  (12) 

o 

Eqs.  (11)  and  (12)  can  easily  be  shown  to  be  equivalent  by  making  axe  of  Eq. 
(6).  Eq.  (11)  cannot  be  integrated  farther  for  x(t)  even  if  one  neglects  the 
R'  term.  Progress  is  possible,  however,  if  we  go  to  the  energy  form.  That 
is,  we  introduce  the  velocity  v  at  a  dependent  variable, 

dx  dv  d2x 

Y  "  dt’  dt  dt7.  (13) 

and  furthermore  use  that 

dv  dv  dx  dv 

-  - - «  V  -  , 

dt  dx  dt  dx 


to  get 


U* 

(L0  +  L'x)v  dv/'ix  +  L'/2  v2  +  2 (R0  +  R_)v  +  2R'  f  xdv  -  b'.  (15) 

y  V 

»  ■*- 

Eq.  (15)  is  In  a  form  suitable  for  a  nontingular  perturbation  expansion  of  the 
solution  where  the  lowest  order  can  be  found  exactly.  For  this  we  introduce 
dimensionless  variables  using  the  principal  parameters  of  the  problem  as 
scales.  The  length  of  the  rails,!,  plays  an  important  role  in  limiting  the 
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final  projectile  velocity.  f#  introduce  a  di vona ionleaa  length  X  by 


-ix. 


(16) 


V#  introduce  a  velocity  acale  Vf  obtained  by  letting  the  initial  aceleration 


.«  2 


L'i 


(17) 


apply  for  the  diatance  £  to  give 

vf  -  V  2aV  -  yT'ii02/.  ' 

Ve  then  introduce  a  dimenaionleaa  velocity  V  by 


(18) 


v  *  VfV  , 

and  a  tine  acale  tf  by 


(19) 


(20) 


A  dimenaionleaa  parameter  h  that  givea  the  ratio  of  rail  inductance  to 
external  inductance  ia 

h  «=  L'i  /L0.  „  (21) 

Another  dimenaionleaa  parameter  r  ia  defined  by 


r  -  2  tf/tc. 


(22) 


where  to  ia  the  time  acale  of  the  external  B-L  circuit, 

Ln 


‘c  " 


(23) 
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Finally  tha  expansion  paraaeter  is  given  by  the  ratio  of  rail  resistanoa  to 
aztarnal  and  plasaa  raaiatanoa 


*o  ♦  *p 


(24) 


In  teras  of  the  diaonsionless  variables  I  and  V  and  the  diasnsionless 
parameters  defined  above,  Eq.  (15)  is  rewritten  as 

V 


(1  +  hi)  V  dV/dX 


+  h/2  V2  +  rV  +  sr  f 


l b*  1  h  /v0tf\ 

where  b  -  — - T  -  -  ♦  -  I— 7“  1 

LoVf  2  2  [  l  J 


XdV  -  b, 

2 


(25) 


votf 

+  r  ■  ;  ■*  -  • 


A  nonsingular  perturbation  expansion  is  set  np  by  snbstitnting  the  series 


V(I)  «  VC(X)  +  eVj (X)  +  s2V2(X)  +  ■ 


into  Eq.  (25)  with  the  initial  conditions 


(26) 


Vo(0)  =  vc/vf  ,Vj(0)  -  0,  i  -  1.2, 


(27) 


The  results  are 


(28) 


(1  ♦  hX)VQ  dV0/dX  +  h/2  V02+  rV0  -  b, 

7° 

(1  +  hX)  (VjdV0/dX  ♦  VQ  dVj/dl)  +  hVoV!  +  rVj  +  r  I  XdV0  -.0,  (29) 

“S  y 

<1  ♦  hlXV.  tV 2 /il  ♦  VidVi/dl  *  VjdVj/dl)  *  h/2  (Vj2  ♦  2V0  V2>  *  rVj  ♦  r  J  XdVj 

O 

for  the  teras  of  order  a°,  a*,  and  s2  respectively. 


0, 
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Propped  Work 

.  -  _  (*o)  is  to  be  solved  analytically  a*  fa*  •• 

Th#  ssqnanoa  of  Eqs.  t2®'  'au/  ” 

».  -1T.I..I  -»•«*-  *•  “  *•  *,,k  ‘ 

solution  of  Eq.  (25). 

H.  ...ly.ic.l  «olot ton  i.  off  to  .  .olid  .t.rt  b.c.n.e  tb.  P.I. 

obtained  .  cto.od  for.  ..Uticn  for  E,  (28)  *■  *»•  »—  1  ‘  I<V»> '  ^  “ 

easily  obtained  by  rewriting  Eq.  (2d)  as 


A  siaple  integration  gives 


^  VodV? _ f  ■  *?-- 

J  b  -(h/2^02-rV~  1  ♦  « 


(31) 


1  +  hi  -  -7 
a 


q2_j  /a  ♦  1  a  -  PO^VilV71 

r-p’/(l+hV0/r)'Z  y  o  ■  i  o  +  p( l+hV0/r) )  , 


(32) 


•bet.  ft  -  (1  *  bvo/tVf)"1  .nd  a  -  (1  ♦  4b fir1)111  . 


(33) 


E,.  (32)  give.  ..  X(V0)  for  .»  id..l  Unneb.r.  Sine.  on.  c.nnot  invert  for 

Vo(r).  ..  .i-plT  ..dlfy  E,.  (29)  to  continne  vitb  the  ..ynptotic 

needed  i.  to  con.id.r  V2,  ete  ..  fnnetion.  of  V0. 


All  that  is 


Budact  i£L  KiPllllPt 

Duration:  Juanuary  1,  19S3  to  December  31,  1983 
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SalsylPJ!,.  fuel  ASjJ  Fyjngg  p IMlilM 

Manual  A.  Huerta,  Principal  Investigator 
1983  -  1984  estimated  salary  135049/9  months 

100%  time  during  1  summer  month  3894 

Graduate  Student  Research  Assistant 

100%  time  during  3  summer  months  2000 

Fringe  Benefits:  17.33%  of  faculty  salary  675 

Total  of  Salaries,  Wages  and  Benefits  6569 


Other  Direct  Costs 
Travel  costs 

to  attend  meeting  of  subject  related  to  this  proposal 

Publication  Costs 

Disc  Drive  (permanent  equipment) 

Software,  and/or  other  equipment  to 
upgrade  microcomputer 


^Indirect  Costs:  60%  of 


250 

200 

550 

136 


mi 


mi 


mi 


tmm. 


•The  Coral  Gables  campus  negotiated  indirect  cost  rate  for  the  University  of 
Miami  for  fiscal  year  6/1/82  to  5/31/83  is  131.0%  of  Modified  Total  Direct 
Costs  (MTDC) .  The  University  is  willing  to  accept  an  indirect  cost  rate  that 
ia  lower  than  the  negotiated  rate.  However,  if  the  proposal  is  awarded,  the 
University  retains  the  right  to  exercise  any  cost  clause  or  cost  regulation 
made  a  part  of  the  signed  agreement. 


